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Abstract 

The vacuum energy of a scalar field in a spherically symmetric background field is 
\Q | considered. It is expressed through the Jost function of the corresponding scattering 

problem. The renormalization is discussed in detail and performed using the uniform 
asymptotic expansion of the Jost function. The method is demonstrated in an simple 
explicit example. 
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1 Introduction 

o . 

The evaluation of quantum corrections to classical solutions play an important role in 
qq | several areas of modern theoretical physics. The classical solutions involved may be 



0\ 



monopoles [|l], §, sphalerons || and electroweak skyrmions j[| ||, [7j 0, ||, [10], [II], [L2|] 



In general the classical fields are inhomogeneous configurations. Thus, as a rule, the 
effective potential approximation to the effective action, where quantum fluctuations are 
integrated out about a constant classical field, is not expected to be adequate. The 
Ph! derivative expansion |13[ improves on this by accounting for spatially varying background 

fields, being a perturbative approximation it has however its own limitations. Having in 
mind that even the classical solutions are often known only numerically, it is clear that it 
is desirable to have a numerical procedure to determine the quantum corrections. Some 
efforts in this direction have already been undertaken [14, 15, [16], 17|. 



X 

The aim of the present article is to develop further a regular analytic approach which 
reduces the evalution of quantum corrections to the corresponding quantum mechanical 
scattering problem. This approach has been developed for theories in (1 + 1) dimensions 



Tq , |i9fl , however, the difficult problem of the summation over the angular momentum 
necessary in (3 + 1) dimensions with spherical symmetry has not been addressed there. 
It is thus the aim of the present article to develop further the regular analytic approach 
and apply it to theories in the physically most interesting (3 + l)-dimensional spacetime. 
Connected with the angular momentum sum, some aspects of the renormalization are also 
more complicated here and are discussed in detail. 
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2 The model and its renormalization 

Let us first describe our concrete model and its renormalization introducing also various 
notations used in the following. We will consider the Lagrangian 

L = !$(□ - M 2 - A$ 2 )$ + -tp(n - m 2 - A'$ 2 V. (1) 

Here the field $ is a classical background field. By means of 

V(x) = A'$ 2 (2) 

it defines the potential in ([!]) for the field tp(x), which should be quantized in the back- 
ground of V(x). The embedding into an external system is necessary in order to guarantee 
the renormalizability of the ground state energy. 
The complete energy 

E[$\=E daa .[$\+E ip [$] (3) 

of the system consists of the classical part and the 1-loop contributions resulting from the 
ground state energy of the quantum field (p in the background of the field $. The classical 
part reads 

E class [$] = ^V g + ^M 2 V! + XV 2 , (4) 

with the definitions V g = f d 3 x(\7&) 2 , V\ = J d 3 x& 2 and V 2 = f <i 3 x$ 4 . Here M 2 and A 
are the bare mass respectively coupling constant which need renormalization as we will 
explain in the following. For the ground state energy one defines [2(J 



Z in) 

where \x is an arbitrary mass parameter and s is a regularization parameter which has 
to be put to zero after renormalization. Furthermore, A( n ) are the eigenvalues of the 
corresponding wave equation 

(-A + V(x))<P (n) (x) = \ 2 {n) <p {n) (x). (6) 

For the moment we assume the space to be a large ball of radius R as an intermediate 
step to have discrete eigenvalues and thus a discrete multiindex (n). 

It is convenient to express the ground state energy @ in terms of the zeta function 

CvOO=£(A 2 n) + m 2 r (7) 

(n) 

of the wave operator with potential V(x) as defined in (|) by 

1 



EM = ztv(s - 1/2)/. 



2.s 



In general £y(— 1/2) will be a divergent quantity and a renormalization procedure for the 
definition of is needed. It is easily discussed in terms of the asymptotic t — > 

heat-kernel expansion associated with the wave equation ([]), 

/ 1 \ 3/2 oo 

#(i) = £ex P (-A^ ~ — e-^^V- (9) 

(n) t-0 V47Ft/ j=0 

The heat-kernel expansion in Eq. (|9|) will also contain boundary terms which in addition 
to the terms written down explicitly will lead to half integer powers in t. However, 
ultimately we are interested in the R — > oo limit and we will subtract the Minkowski 



r 



-2- 



space contribution in order to normalize E v \f^ = 0] = 0. For the case that V(r) 
for r — > oo, e > 0, these terms will then disappear. 

For the renormalization, only the first three terms are relevant which read explicitly 
A = J d 3 x, A 1 = — J d s xV(x) and A 2 = (1/2) / d 3 xV 2 (x). Their contributions to the 
ultraviolet divergencies in the groundstate energy read 



1 64tt 2 \s m 2 2 ' 

+^f 1 + '»^-iK (10) 

1 ^ ■ i,/" 2 
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327r 2 ys m 

The first term is a constant independent of the background field $ and it can be dropped 
(especially it will be absent after subtraction of the Minkowski space contribution). In a 
more general context of a gravitational background field this term would yield a renor- 
malization of the cosmological constant. The second term can be absorbed in a renormal- 
ization of the mass M of the background field 

9 9 X'm 2 ( 1 , m 2 \ , . 

- M + i^(-; + 1 + ln vj 

and the third one in the coupling constant A by 

x A' 2 / 1 , m 2 \ 

647T 2 \ s 4/j 2 J v 7 

The kinetic term V g in E[$>] suffers no renormalization. By defining 

E%" = E V [*\-E* U (13) 

one obtains the finite groundstate energy, which is normalized in a way that the functional 
dependence on $ 2 present in the classical energy is now absent in the quantum corrections 

rpren 

E v • 

Let us note that this is just the well known general renormalization scheme written 
down here explicitly in the notations needed in our case. 
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3 Scattering theory and Jost functions 



Let us now restrict to a spherically symmetric background field $(r). Then the multiindex 
(n) — > n, I, m consists of the main quantum number n, the angular momentum number / 
and the magnetic quantum number m. In polar coordinates the ansatz for a solution of 
the wave equation (§) reads 



= -<l>n,i(r)Yi m (6,(p) 



where the radial wave equation takes the form 

" d 2 i{i + r 



dr 2 



V{r) + K, 



0n,z(r) = 0. 



1; 



(2) 



Now we use the standard scattering theory within r G [0, oo) and have the momentum p 
instead of the discrete \ n j. Let 4> P) i{r) be the so called regular solution which is defined 
as to have the same behaviour at r — > as the solution without potential 



Jiipr) 



(3) 



with the spherical Bessel function ji |?T). This regular solution defines the Jost function 
fi through its asymptotics as r — > oo, 



i r 



fi(p)h t (pr) - fiip)hf{pr) 



(4) 



where h^(pr) and hf[pr) are the Riccati-Hankel functions 

Now we use the Jost function to transform the frequency sum in Eq. (|7|) in a contour 
integral. Let us assume for a moment that the support of the potential is contained in the 
cavity of radius R. Then the above Eq. (^) gets exact at r = R and may be interpreted as 
an implicit equation for the eigenvalues p = \ n ^. Choosing Dirichlet boundary conditions 
at r = R, 4>p,i{R) = 0, it reads explicitly, 



fi(p)hr(pR)-K(p)h+(pR)=0. 



(5) 



As already mentioned, ultimately we are interested in the limit R —>■ oo and in that 
limit the results will not depend on the boundary condition chosen, once we assume that 
V(r) ~ r~ 2 ~ e for r — > oo. 

Let us now consider the ground state energy associated with the eigenvalues deter- 
mined by (H). It is convenient to represent the frequency sum in @ by a contour integral, 
the basic idea being explained in detail in for example pZ| , [Z5| |. Using Eq. (|5|), one imme- 
diately finds 



EM = M 2s £(/ + l/2) [^(P 2 + 

oo 

V s E(^ + V2)E(^ 2 -« 



m 



2\l/2-s 



^-HMpfhripR) - fttp)ht(pR)\ 



n.lj 



.1/2- 



(6) 



1=0 
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with — k 2 i as the energy eigenvalues of the bound states with given orbital momentum 
I. The contour 7 is chosen counterclockwise enclosing all real solutions of Eq. ([5]) on the 
positive real axis. The division of the discrete eigenvalues within the large ball into positive 
(inside 7) and negative ones (— is determined by the conditions that V(r) — > for 
r — ► 00. In that way, in the limit of the infinite space the negative eigenvalues become 
the usual boundstates and the X n> i > turn into the scattering states. 

For the calculation of (H), as the next step one deforms the contour 7 to the imaginary 
axis. A contour coming from ioo + e, crossing the imaginary axis at some positive value 
smaller than the smallest K n and going to ioo — e results first. Shifting the contour over 
the bound state values K n which are the zeroes of the Jost function on the imaginary 
axis, the bound state contributions in Eq. @ are cancelled and in the limit R — > 00, 
subtracting the Minkowki space contributions, one finds 

E V [<S>} = -^// 2s £(Z + l/2) J dk [k 2 -m 2 ]^ s ^lQ/,(i*). (7) 

i=o m 

This is the representation of the ground state energy (and by means of (H) as well of the 
zeta function) in terms of the Jost function, which is the starting point of our following 
analysis. It has the nice property, that the dependence on the bound states is not present 
explicitly, being however contained in the Jost function by its properties on the positive 
imaginary axis. To the authors knowledge, this representation of the ground state energy 
respectively for the zeta function is not known in the literature. Its analog proved to be 
useful for explicit numerical calculations already in the case of a potential depending only 



on one coordinate ||19|| . It is connected with the more conventional representations by 
means of the analytic properties of the Jost function. Expressing them by the dispersion 
relation 

«*)-n(i-^)-p(-|r^(«)). 

where Si(q) is the scattering phase, we obtain from ([/]) 




-l-2s 



m — \ m' — k 



dq , 2s 5i(q) , 

\/q z + m z ) 

which gives the expression of the ground state energy through the scattering phase. From 
here one can pass to the representation through the mode density by integrating by part. 
Note, that this representation can be obtained also directly from (^]) in the limit R — > 00 
by deforming the contour 7. 

Let us add a discussion on the sign of the ground state energy. In (|5|) the first contri- 
bution results from the bound states and is completely negative. The second contribution 
which containes the scattering phase Si(q) is positive (negative) for an attractive (repul- 
sive) potential, i.e., for V(r) < (V(r) > 0) for all r 0. So the regularized (still not 

3 This is a well known fact from potential scattering, see e.g. 
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renormalized) ground state energy (BI) is positive for a potential which is repulsive 

for all r (there are no bound states in this case) and it is negative for a potential which is 
attractive for all r. Now, if we perform the renormalization in accordance with ( p]) and 
(H^), we obtain 

E ren = E v [$\ + ^ ( - + log ^ - 1 ] j°° dr r 2 V(r) 
8n \s m 2 J Jo 




drr 2 V(r) 2 . (9) 

The contribution containing A had been already subtracted in £^[3>] (0). This expression 
is finite for s — > 0, i.e., when removing the regularization. But due to the subtracted terms 
there is no longer any definite result on the sign. Note, that this is in contrast to the case 
of a one-dimensional potential where it had been possible to express the subtracted terms 
through the scattering phase [19| . 



4 Uniform asymptotics of the Jost function 

Let us continue with a detailed analysis of the ground state energy, Eq. (]?]). As we 
have easily seen using heat-kernel techniques the non-renormalized vacuum energy 
contains divergencies in s = 0, see Eq. (|T0|), which are removed by the renormalization 
prescription given in Eqs. (|ll|) - (0), resp. (|). The poles present are by no means 
obvious in the representation (|7|) of £^[3>]- However, in order to actually perform the 



renormalization, Eq. (|13l) , it is necessary to represent the groundstate energy Eq. (0) in 
a form which makes the explicit subtraction of the divergencies visible. This will be our 
first task. 

As is known from general zeta function theory as well as one sees from simply counting 
the large momentum behaviour, the representation Eq. (|7]) of i^[3>] will be convergent 
for Res > 2. However, for the calculation of the ground state energy we need the value 
of Eq. (^|) in s = 0, thus an analytical continuation to the left has to be constructed. The 



basic idea is the same as the one presented in p2| |23|] : adding and subtracting the leading 



uniform asymptotics of the integrand in Eq. (0). Let 

E v [$] = E f + E as (10) 

where 

COS(TTS) 2s ^ °f „2li-« 9 



E f = -2^2^^ + 1/2) J dk [k 2 -mf2-^[l nfl (ik) -In fr m (ik)\ (11) 
i=o m 

and 

^ = -^^^ 2s f:( ; + i /2) Jdk [k 2 ~m 2 }^^-\nfr m m. ^ 

7r J ok 
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The idea is, that as many asymptotic terms are subtracted as to allow to put s = in the 
integrand of Ef. This term will then (in general) be evaluated numerically. In E as the 
analytic continuation to s = can be done explicitly showing that the pole contributions 
cancel when subtracting Eq. (fLOf). Note, that the contribution resulting from Aq 

has been dropped already in Eq. (0). 

The first task thus is to obtain the asymptotics of the Jost functions. This may be done 
by using the integral equation (Lippmann-Schwinger equation) known from scattering 



theory [21]. For the Jost function one has (y = I + 1/2) 

oo 

fl(ik) = 1 + Jdrr V(r)^ lk (r)K u {kr), (13) 
o 

with the regular solution given by the integral equation 

r 

<PiAr) = h(kr) + / dr' r' [I v {kr)K v {kr') - I u {kr')K u {kr)]V{r')<p lAk {r'). (14) 



o 



General zeta function theory tells us that the divergence at s = contains at most terms 
of order V 2 . Thus one might expand In fi(ik) in powers of V and take into account only 
the asymptotics of terms up to 0(V 2 ). The expansion in powers of V is easily obtained. 
Using Eqs. (|IB|) and ( |T4"D one finds 

oo 

In fi(ik) = dr rV(r)K v (kr)I v {kr) 



o 

oo r 

dr rV(r)Kl(kr) J dr' r'V(r')I 2 (kr') 
o 

+0(V 3 ). (15) 

Now the uniform asymptotics for / — > oo of In fi(ik) is essentially reduced to the well 
known uniform asymptotics of the modified Bessel functions K v and I u , (9.7.7) and (9.7.8) 
in H. With the notation t = l/yjl + (kr/u) 2 and rj(k) = yfl + [krjv) 2 + \n[{kr/v)/{\ + 
1 + (kr/u) 2 )], one finds for v — > oo, k — > oo with kjv fixed, 

I v (kr)K v (kr) ~ -L + ^fi-Gt 2 + 5^ + 0(1/^) 

WK , (ir) ~ ^ (1 + (tr/y)ami + (tr , /ym J i + 0(iA)]- 

Using these terms in the rhs of Eq. fllBD we define 



1 °° 

\nff sym (ik) = J dr — 



r V(r) 



n * 1/2 



7 



8v 3 



dr 



r V (r) 



3/2 • 



(16) 



Thereby the r'-integration in the term quadratic in V has been performed by the saddle- 
point method using the monotony of f](k). Now by means of flH)| ) the limit s — > can be 
performed in Eq. (|TTJ) and we obtain 



E f = ~ E( ; + 1/2) / dk v 7 ^ 3 ^^ Q*fi(ik) ~ ln/r m (^)) , (17) 

a form which is suited for a numerical evaluation. 

For E as at s = one might explicitly find the analytical continuation. First of all, the 
/c-integrals may be done using 



dk \k 2 -m 2 }^- 8 -^- 
ok 



kr 



i+Kf 

\mr J 



r(n/2) 



to yield Eq. ( |A.1| ) in the appendix A. ^From this, the renormalization, Eq. ( |T3"D can be 
carried out and we arrive at 



TP — Tpd™ - ___ 
v 487T 



drV{r) 6mV 2 - 12Cr(-1) - 31n2 



6r 2 ln(16mV) m 2 + -V(r) - 7(1 + 12m 2 r 2 + 6rV(r)) 



+ E, 



sum 
as 



The complete energy E[$\, (§), consisting of E dass [^}, (|), of the background field $ and 
of the renormalized loop contribution (|i~3"D 



rpren 



E^]-E d ™ = E f + E as -E" 



reads 



E[*) 



1 V g + \m 2 V x + X ren V 2 



2 y 2 

00 



(19) 



(20) 



48tt 



drV(r) 6m 2 r 2 - 12C^(-1) - 31n2 



-6r 2 ln(16m 2 r 2 ) ( m 2 + -V(r) ) - 7(1 + 12mV + 6r 2 \/(r)) 
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The remaining task for the analysis of the ground state energy in the presence of a 
spherically symmetric potential is the numerical analysis of the above quantity. To achieve 
that, a (as a rule numerical) knowledge of the Jost function fi(ik) is necessary Apart 
from this only integrals of the potentials and convergent sums have to be dealt with which 
presents no problem. 



5 Example: Square-well potential 



If the potential has a compact support, there is a formalism how to obtain the Jost 
function (at least in principle). Starting point is the observation, that one may write the 
regular solution in the form 

Mr) = ui, p (r)e(R -r)+ % - [fi(p)K(pr) - f?(p)ht(pr)} 6(r - R), (21) 
where V(r) = for r > R is assumed. The matching conditions then read 

MR) = \[fiip)kipR)-f:{v)ht{vR)}, 
M R ) = ^p\fi(p)K'(pR)-fttp)ht'(pR) ■ 



Combining the two equations and using the fact that the Wronskian determinant of hf 
is 2i, one arrives at 



hip) = -~ (p Ul , p (R)ht'(pR) - u' ltP (R)ht(pR) 



(22) 



which gives the expression of the Jost function for a potential with a compact support 
through the wave function. 

For the square- well potential V(r) = VqQ (R — r) it is easily seen that 



ui p [r) 



p 



i+i 



with ji the Riccati-Bessel function and q = y/p 2 — Vq. 
So we obtain the well known formula 

friik) = R (j) [qI' v {qR)K v {kR) - H u (qR)K' u (kR)], 



(23) 



with q = yjk 2 + Vq. This has to be used for the numerical evaluation of Eq. Q2"0|). Due to 
the simple form of the potential V(r) the r-integrals may be done explicitly and numeri- 
cally easy tractable expressions result. For Ef we obtain 



1 00 7 
E f = V(/ + 1/2) / dk [k 2 

71 2-0 J 



(24) 
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in frm 



2k 2 



V 



-1/2 



V 2 



8uk 2 



R 2 + 2{vjk) 2 
(1 + (kR/v) 2 ) l/2 " VA; 



- 2 



21? 



1/2 \ 
- 1 



V 



-3/2 



+ 



1 + 



-5/2' 



which is, as one might easily check, a finite expression. For the renormalized contributions 
of the asymptotic terms, we obtain 



rp _ rpdiv 

11 'as J^tr, — AO 

* 487T 



V R 



-12Cr(-1) + ^-m 2 R 2 -3 In 2- 2m 2 R 2 \n{lQm 2 R 2 ) (25) 



+-V R 2 - V R 2 ln{16m 2 R 2 ) - 7 (1 + 4m 2 R 2 + 2V R 2 ) 
3 



I rpsum 
' as,swi 



where E s ^ w is given in appendix B. It is presented in form of quickly convergent series 
which allow an easy numerical analysis. 



E r v en l®} V =0.9 V Q = 0A Vb = 0.1 




V = -0.4 

Fig. 1 The complete energy as a function of the radius R for different values of the 
height of the potential wall V for m = 1. 
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The results for E™ n [<&\ are presented in Fig. 1. For small values of the radius R of the 
support of the potential, it is seen that even for negative values of the potential, V < 0, 
a positive vacuum energy results. The reason is, that for small values of R the bound 
state energies are located closely to zero, giving only small negative contributions. For 
increasing R, their number and their values are increasing leading necessarily to negative 
vacuum energies. 

The behaviour described is the one we expected from the one-dimensional results 



presented in [19]. However, the absolute orders of magnitude are by two lower than in 



the corresponding one- dimensional considerations. 



6 Conclusions 

In this article we reduced the task of calculating the vacuum energy of a scalar field in 
the presence of a spherically symmetric background field to the corresponding quantum 
mechanical scattering problem. We were able to present the renormalized vacuum en- 
ergy solely in terms of the quantum mechanical scattering dates summarized by the Jost 
function. For the example of the square-well potential we showed that a direct numerical 
analysis of the vacuum energy is possible. 

Several extensions of our approach are necessary. From the physical point of view 
the consideration of higher-spin fields is necessary and envisaged. In addition, in order 
to apply our formalism to classical background fields like f.e. sphalerons and electroweak 
skyrmions, a numerical analysis of Eq. (O) is necessary for cases when the Jost function is 
known only numerically. However, the scattering theory developed during the last decades 
provides many techniques and results so that also here progress seems possible. 
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A Contributions of the asymptotic terms to the zeta 
function 



Here we calculate the analytic continuation of E as (12) using the asymptotics of the Jost 
functions \afi s {ik) After carrying out the k- integration by means of (pl|) we obtain 



E* 



T(s)ii 



™ oo 
°9 oo 

/ dr rV(r) E 



v 



\mr 



m 1 2s fj? s 



2^T(s - 1/2) 



r(a + i) . 

q I 4(mr) 3 



r(» + 2) tJ 



mr 



\mr J 



-s-2 



3(mr) 



\mrj 



-s-l 



-l -s-3 



\mr J 



2s 



4v^Fm 2s + 2 r(s - 1/2) 



/ dr K 2 (r) ^ 



i + i- 

mr 



l s-l 



(A.l) 



The analytical continuations of the above expressions to the value s = is easily obtained 
and listed below, 



E n 



— I dr V(r)(l - 12mV) (- + In r 2 /i : 
96-7T i V s 

dr V(r) \-l - 12Ca(-1) + 12m 2 r 2 (l - 7 -31n2) 



48tt 



oo _^ 1 °° 

g — y dr V(r) f- + lnrV) + 4^/ rfr ^0)[1 - 7 - 31n2] 



oo 1 °° 

— J dr r 2 V 2 (r) (- + In r 2 /i 2 ) + ^ J dr rV 2 (r)(l - 7 - 3 In 2) 



167T 



(A.2) 



where the sums in 



Tpsum 



I r 

- dr ^(r)E 
1 ? 

+-/dr 7(r)E 



/ mr 
\ v 



M 111 
v I 4 



-In 1 + 



f mr 
1+1 — 



1+ 



\ P J 



2 

+ 3 



1+1 — 



1 

12 



mrV 



12 



i 7 00 l 

— - / drr 2 V\r)Y- 

8vr J ,n" 
n f=0 



fmr\ 



- 1 



(A.3) 



converge. 

This expressions have been the basis to give the results for the special example listed 
in appendix B. 

B Asymptotic contributions for the square well po- 
tential 



In this appendix we give the result for E™™ w which reads 



rpsum 
as,sw 



V R™ 

> V 

Atx ^ 

* n 1=0 

V n R°° 



1 /mi? 
3 V V 



-In 1 + 



mR 



+ 2- 



2u 
Rm 



arctan 



Rm\ 



v 



47r ir ^ 

V*R 3 ~ 1 



1 / (mR 

i 1+ — 



mi? 



mi? 



6 



arctan 



1 ( (mR\ 



1 

12 
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These sums have been used for the numerical analysis of the groundstate energy in the 
presence of a square-well potential 
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